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The work consists of tuo parts . In part ’A’ a 
parnnetr ic study has boon nade on the condensation polymerizatioi 
of A3 LB type '.’.or. oners and on polymerization between monomers of 
typo A. .A and. Ail’l, for the caso where tlie equal reactivit]- 
hypothesis is rot valid. The analysis has been done for a 
segregated OSTk. For Af.B type monomer the molecular weight 
distributions and their moments have been compared for two 
different kinetic models, suggested before. These bwo models 
have been ext ended to the other case and a third model is 


proposed. 

In part ‘B 1 reversible hydrolytic polymerization 
of e-caprolactam to nylon 6 in a batch reactor has been 
studied. The primary purpose was to simulate the reactions 
in the polymerization process and to establish a computationally 
efficient mathematical technique for solution. The approach 
was not successful, however, for the system considered here 
due to a problem of computational unstability. Possible 
reasons for the failure and possible ways to overcome this 
problem have been discussed. 
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Part A. 
fPO, fQO 

k I 

K 

p 

M 
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fractions of initial monomers P^nd , 
respectively in ARA + BR'B 
the first rate constant in the models of 
unequal reactivity 

the second rate constant in the models of 
unequal reactivity 

total original concentration of monomers in 
ARA + BR'B 


dimensionless k-th moment in ARB type condensation 
polymerization 

=(k) =(k) =(k) 

M p * ^-th moments for distributions of species P,Q & S, 

respectively 

P 10 , ^10 ,S 1 0’ Initial monomer concentrations in ARA 4- BR'B 

polymerization. 

? n bifunctional chains with both end groups A in 

ARA +■ BR'B type polymerization 

Q n bifunctional chains with both end groups B in 

ARA + BR'B type polymerization 

R ratio k-r/k 

■*“ ir 

S n bifunctional chains with A & B as the end groups, 

in either ARB type of ARA + BR'B type polymer izat it 
dimensionless batch time, k^P^ q t for ARB, and 
k p Mq t for ARA + BR'B 


X 



n dimensionless concentrations in ARB type 

polymerization 

^ actual reactor - holding time in OSBR 

® dimensionless reactor holding time in CSTR, 

k p P 1q T for APB and k for ARA+BR'B 

^IC f / l -u Number average and weight average chain lengths 

Polydispersity Index 
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B 
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n 

G(x) 


G 3 

=k 

G 

x 
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.° v° 

u o * JEV-rr 


k c k c k c 


initial concentration of acids or amines 
concentration of monofunctional chains of length n 
constant used in the~ eqnatrohhdepleting' 1 the 
functionality of cyclization equilibrium constant 
K 3 (n) with n 
monomer concentration 

concentration of macrorings with n repeating units 
coefficients of laguerre polynomials 
the normalized distribution 
abbreviations for the j-th moment of G(x) 

k-th moment of normalized distribution G(x) 
rate constants for the uncatalyzed reactions in 
caprolactam polymerization 

rate constants for the catalyzed reactions in 
caprolactam polymerisation. 



‘^2 , ^ c ' 


k- ,k. 



K r K ? 


K 3 (n) 


L (x) 
n 

M(v) 


= (k) 
M 

V 


v 

V 

X 


r (x) 
0. 


forward reaction rate constants 

reverse rate constants for ring opening and 

polycondensation, respectively. 

cyclization reaction rate constants as function 

of ringnlength n 

Equilibrium constants for ring opening and 
polycondensation, respectively 

cyclization equilibrium constants, as function of 
ring length n 

generalized la guerre polynomial of order n 

the contineous distribution used for the actual 

discrete distribution 

k-th moment of distribution M(v) 

the number of moments used in Min’s method 

concentration of amino cagroic acid 

concentration of the bifunctional chains of 

length n 

the argument in the contineous distribution M(v) 
water concentration during polymerization 
the argument in the normalized distribution 
function G(x) 

gamma function of argument 'x' 

free parameter in the generalized laguerre 

polynomials 




n-th moment for monofunctional chains distribution 
n-th moment for bifunctional chains distribution 
n-th moment for macroring distribution 
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CHAPTEB A-1 


Int roduction * 

Usually in the study of the kinetics and the molecular 
weight distribution of linear condensation polymerizations, 
the assumption is made that the reactivity of a functional 
group * is independent of the nature and length of the chains 

( i _5 ) 

to which it is attached . 

However, many systems have been reported for which this 
assumption is not valid, e,g. polymerization of alkali metal 
salts of p-chloro- and p-fluorothiophenols , formation of 
polyethylene terenhthalate^ ^ and polymerization of 
diisocyanates^’ 10 ^ and of phenol for maldehyde systems 
The violation of the equal reactivity assumption for polymeri- 
zation between acid anhydrides and aromatic diamines has been 
reported ^ ^ ^ , and it has been shown^^ that the MWD in the 
condensation product of divinyl benzene with p-cresol is 
different from that predicted using this hypothesis. 

Model studies on this change of reactivity have been 
relatively fewer. Kronstadt et.al/ 16 ^ used two different 
rate constants in their kinetic scheme with partial success 
in matching the theoretical prediction with experimental data. 
Miller and Macosco^ 1 ^ have considered stepwise polyfunction 
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polymerisation where functional group reactivities in the 
same monomer change due to a substitution effect. 0? Izmir and 
Odian have analysed linear step polymerization between AA and 

(iq) (20) 

tC monomers'' y) and between AA,BB and CC type monomers , 
where B and C are functional groups of the same nature but 
having different reactivity. Gandhi and Babu have considered 
the change in reactivity of functional group due to asymmetry 
or substitution effect in the molecule^ 2 ^, and due to origi- 
nally present asymmetry and induced asymmetry on the same 
molecule . Kumar, Gupta and coworkers have considered 
violation of the equal reactivity hypothesis in phenol 
formaldehyde systems^ •'’-15 ) an £ ^ ave proposed two kinetic 

( 21 - 23 ) 

models on general linear polycondensation of ABB type monomers 
However, all these are studies on batch reactors. , 

Condensation polymerizations in CSTB's are important 

, ( 26 ) 

industrially but few researchers have worked on it. Biesen teerger 

computed the polvdispersity indices in a micromixed model of 

„ ^ ..( 27 , 28 ) 

CSTB . Based on the same hydrodynamic model, Gupta et ai 

obtained the MV-D and its moments for the ..-. < two proposed kinetic 

( 21 - 23 ) 

models by them earlier for batch reactors 

It should be noted that for the polymerization in a 
CSTR, a. microscopic segregation model .is often closer to the 
actual case than the micromixed model particularly because of 
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( 29 ) 

the high viscosity involved. Tadmore and BiesenHierger 
found the MWD’s and their moments for condensation polymeri- 
zations carried out in a segregated CSTR, for ARB type monomer 
fo 1 lowing the equal reactivity hypothesis. In this work we take 

/ /p A } 

the two models proposed by Gupta et al ' and make a 

theoretical parametric study on the corresponding MV/D ' s and 
moments in a segregated CSTR. We then extend these models which 
were originally for ABB type, to ARA + AR'B type polycondensations 
and also pronose a third model to make a closer approximation 
of some real cases. In both the cases the reaction is assumed 
to be irriversible which implies that the condensation product 
is withdrawn as soon as it is produced. 
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CHAPTER 1-2 
PCPMITLATIOF 

AR Bjty pe monomer : In such a reaction at any time only one 

family of molecules 'will be present, namely, S . The general 

n K n 

reaction is: S m + S n — S m+n + condensation product. The 
two models proposed are the two limiting cases of real 
polymerizations : 



The corresponding rate equations for species and moments, in 
dimensionless form, for batch reactors, are given in Table A1 . 


In an SCSIR there is a residence time distribution and 
in an ideal CSTP, it varies from zero to infinity with a RTD 
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function E(t) dt = i exp(-t/'t‘ )dt 


(30) 


(3) 


(29 30) 

The average outlet concentration of S n is given by ’ 
oc 

S n i ijs n (%) exp(-t/ / T )dt (4) 

O 

To make the equations dimensionless we define x = k^S^t a:. 

arid 6=kS 10 t , (5) 


S 10 


So the d imensi onless outlet concentration from the segregated 
CSTR is >"iven by 


A 

s 

2 _ __ _jl 

n S 


10 


1_ 

e 


r 


z n exp(-x/£)dx, n = 1,2,... (6) 


vJ 

Q 


Prom the batch concentrations of S n one can get the outlet 
concentration by integrating equation (6). Molecular weight 
distributions may be calculated by the following equation 


mass fraction of n-mer 


/v 

nS 


n 


A 

nZ 


InS 


n 


( 7 ) 


n 


The values of the moments in the out let of the SCSTR can 
similarly obtained by using the equation 


A 

"k 


= ~ eXP k=0,1 , 2 ... 


where n ^ ^1 0 » k= 0,1,2 


( 8 ) 

(9) 


The molecular weight averages can then be obtained as follows 


K = % / s c> Y w= V*s > f = /van 


( 10 ) 



b 


ARA + BR ' B type monomers i At any time in the reaction mass, 

three types of species >1111 he present viz. P n , the odd A-mers; 

Q . the odd B-mers; and S , the even-mers, as identified originally 
n n 

by Plory^\ Bach species will have a distribution of its own. 

The moments for the distributions are defined as follows 


= (k) 
M 

P 


£ (2n-l) k P *, 2 (2n-D k Q 5 

11 y - - 


=(k) 


=(k) 

M 






I <2»> 

* 0*1 


k 


n 


The reaction scheme for the three models are shown in 
Table A-2. The rate equations for various species and moments are 
shown in Tables A-3.1, A-3.2 and A-3.3 for batch reactors. The 
values at the reactor outlet will be obtained as in equations 
(6) and (8) above, the only difference is that f x f and ' ’ now 
defined as x = k p M q t & 9 = k p M Q ^r ....(12) 

In Model I, we assume that the monomer-monomer reaction 
rate constant is different from those for all other reactions, 
i.e. between monomers -polymers and polymers -polymers . In Model II 
we are assuming that the reactivity change due to the substitution 
effect is same for both the functional groups. In Model HI we are 
assuming that the reactivity of one of the functional groups 
changes due to the substitution effect, but that of the other 


remains unaffected. 
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Table A-l 

Rate equations in ARB type polymerization 


Model I 


Model II 


dZ,/dX = -(R-1)Z“ - M Z- 
I 1 o 1 

d V dx = § zf - m o z 2 


dZ x /dX = - RZ-jMq 


dz 2 /dx - 2 Z i (R“1)Z 1 Z 2 - M q ^ 2 


n~l 


dz/dx = x y~ Z Z 
n 2 f—; m i 

m=l 


n -1 

- M Z , dZ /dX = (R-l)Z-Z - + ~ Z Z 
n-ra o n' tv 1 n~l 2 m 


n 


dM 

o 

5 x~ 


n — 3 g A: f m « * 


(R-l) 7 2 M o 

2 2 


dM 

o 

a>r 


*“(R— 1 )Z- Z ~ M Z , 
In on' 


n = 3, 4, . . . 

¥ z ? - 5 s - 


dM k 

dX 


(2 lt " 1 -l)Z^ 


1 —IS 1 

+ 2 lTTk^irr M l M k-3 


- M k M Q/ k = 1,2, ... 


dX 


(l»2 k “ 1 )(R-l )zj + (R-DZj 


k-1 


^ li (k-»l )i M 1 
1 k‘ 

\ £ nrfenr M i M k-i 


Initial conditions ; at X = 0: Z^ = 1, Z 2 = Z 3 = * * * Z n = °* 

M^. — 1 f k = 0, 1 , 2 . . . 


Note: For definitions of the symbols see Nomenclature. 



Table A-2 


Kinetic scheme for various AHA + BR'B polymerization 

k q 

General reactions: P + Q -EliBSs . > + ¥, n,m = 1,2,... 

n m ^n+m-1 * * ’ 

1c 

P + s P L + ¥ , n,m = 1,2,... 

n in ^ n+m 9 9 

k 

Q + S : n >- Q , + W , n,m = 1,2,... 

m ^n+m ’ * 


n 


k 


S + S '. s . a ±£^ S ^ +¥ , n,m = 1,2,.. 

n m ^ n+m ’ ’ 


Model I 


k pQ 

,H “ 

k i 

; k pQ 

, In 


S’ 

n - 2,3, ... ; k pQ>n1 - 




k r\ 


, 

k ; 

n, m = 2 , 3 * 




pQ 

,nm 


P 


k„ 0 

ps 

,nm 

k p- 

n,m 

= 1, 

, 2 , 

9 * • • 


k Qs 

,nm 

k p’ 

n,m 

= 1, 

>2, 

M » • 


k ss 

,nm 

k p’ 

n ^ 

m,n, 

,m 

= 1 

, 2 i . • » 


= k p /2, n =m, n,m = 1 ,2, 


Model II 


Mod el III 
v 


k 

k 

k 

k 

k 

k 

k 

k 

k 


pQ,1m 


k i ; 

m = 

1 , 2 , « • * 5 

k pQ ,n1 k l* 

n= 1,2,... 

pQ,mn 

= 

k p ; 

m,n 

= 2,3,*. 

. 


ps ,1m 

= 

k i ; 

m = 

1,2 , ♦ * • ; 

k ps,nm = k p’ 

n=2,3, . . ;m= 1,2,. 

Qs ,1m 

= 


m = 

1,2,..*; 

k Qs ,nm = k p ; 

n= 2,3,..» m= 1 , 2 , 

'ss ,nm 

= 

k: 

P 

n ^ 

m, n,m = 

1,2,*.. 



= 

k p /2 , n = m, n, 

m =: 1,2,..* 



pQ,1m 

= k I> 

m= 1,2,...; 

k pQ,nm k 

ps , 1 m 

*s 

M 

X 

II 

m= 1,2,...; 

k = k. 

ps,nm 

Qs ,nm 

II 

* • 

V. 

3 

II 

■>* 

f\> 

V. 

* 

• . 

ss ,nm 

■ V 

n 4 m, n,m 

= 1,2,... 


= k ,n=2,3, ... ; m= 1 , 2 , 



Table A3.1 
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Rate equations for Model I in ARA + BR* B polymerization 


i = - k ( R - l ) P lQl . kF ^' 1 + 


dP 
dt 

dQ, . . 

ar = - k P ' R - 1 ’ p i°i - * p °i <® p + “ s ’) 

dP n "*^- 

®T = k p Z P Ar - V „< +g “> » 11 = 2,3,... 


P /— _ r n-r "p" 

r=l 


dQ 


n 


n 


Ht“ = k p 


H P r Q n-r+l + ¥ ^p 


n— 1 

k p E s - s 


r=l 


r=l 


r ^n— r “ If k p S n 


s(o> 


kpSp < H < o) +^< o) ) ', 


n — * 2 ^ 3 / m * m 


aM ^ o) 

dt 


dffi 


(o) 


at - = -*p (R - 1J p i Q i - k P ^ o) ^ o) > 


= kp (R-l)P 1 Ql + k p M<°> t^ o) - i k p k p H<°> 


( H ^ o) ^ 0l ) 


d ft 


(k) 


-k p (P-DPi Q X - k p H< k) Hq° t k p £ 1 ;(Li)t ff p 1I ~ 1) g, 


1 = 1 


Ft — 

at 3 k p ( 8 -i) p iQi -k p ^ k) k p £ XTTCTT ! 

d $ i k) 


kl jBr(k-l) ff (l) 

S 


_k_ 

L 

1=1 


kZ _^( k - l ) M ( l ) 
Q S 


FF 


2- kp(B-i) PjOj. + kp £ TiTferTT “q 


Ef<l) g(k~l) 


1=0 


. 1 kl S {1) sf(k— 1) 

+ 1 k p /- rrn ^ rrr ” R ~ K ~ 

1=0 


-4 k p 0 j o) - k p H < k) < ff ' o) + fi £°>) 


Continued . . . , 



Table A3«l (continued ) 


at t = 

ll 

t-4 

•* 

o 

n 

P 10 » Q 1 

_ Q 1q ; S 1 

P-/ Q_; S = 0 
n n' n 

/ n 

85 2 g 3/ « * * « 

II 

. — s 

o 

'Ifef* 

<4 1) = i. 

. =^ k) = 

:P 10> 

M^°^ = 
Q 

Q " " 

•= = 

Q io ? 

«i 0) : 

= , . 

• = Pii K> 

= o* 


= Oj 


Table A3 ,2 

Rate Equations for Model II 


err- = - R p i ^ o) - R p i 3 i o) 

sr = - R Q i ^ o) - R Q i fl s 0> 


dP 


n— 1 


XT = k p (R " 1) P l S n-l + k p Z P r S n-r ~ k p (lul)P n Q 1 ~ k p P n 


P 

•k P_ 
p n s 




(o) 


dQ 


r=l 


n-1 


err- - k p °1 S n— 1 + k p 2 Q r S n-r - k p ^ Q n P 1 - K p Q n 

(fl^ + fi< o) ) 

o ) - __ 

= k p (R-l) PjQj, - k p (R-l) - k p (R-l) P 1 fi^ o) 

-^^ 0> ^ 0> 

£r~ = k p (R-l) P lQl - k p (R-l) S^\ -kp (R-l) Q^ 

- *p ^ o) s' o) 


(o) 


Continued „ # . , . 



Table A3«2 (continued): 
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dtt 


(o) 


HT - » k p (a-l) P x f^ o) - k p (r- 1 ) p 1 Q x + k_ (B-i) B< o) q. 


+ k 1 


yP 1c _L_ / i'ij-, JT - j + _ 

= P O P OX P S 


( o ) ' 


- k ( R— 1 ) H, 


(o) T 


k 


ct(°) . s(o)' 




^0° _|- k r, S, 


l o) °i 


Ht 


(k) 


r(l) 


- kptR-i) p lQl + k p (R-i) Pj irferr "<r + 

1=1 


k 

k p Z 


1=1 


TK^TTT ®i 1> ^ k " 1) " k P (R - 1> «p k>Q l - k p ^ p i S q 0) 
- k ftl k) rai o) 


-P <¥ “Q 


d^ k) 

31 


= k p dui^ Qj. + k p (R-i) q x x 1 ., k : 1) , ' 


Ft 


( 1 ) 


1=1 


+ k P £ Ti - rfc-T TI s i 1J ^ k_1> - k p <^ l) P 1 s Q k, - k , 


i3r( O ) 


CR— 1) QjH£ u/ - k p H. 


(o) B (k) 
Q 


*6 


(k) 


3F 


• 2 \ t 8 - 11 p i°i + k p P 1 Z rrffciTl ^ 1>+ k 


k_ 

r 

1=0 


ijrztn )m q 1) + k p (IU1) Q i 51 u(k-Di 

• L=u 1=0 

+ * k p £ TTTfcr7r M i 


p 

(1) 


T ^-rrr ^ X) ^ k ~ 1) - | k p F^ k) H^ o) - k p (R~l) 


*£ k> (P 1 + Q l ) - k P ^p 0) + r q 0)) 


Initial conditions : Same as in Table A3.1, 



Table h 3.»3 
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Rate Equations for Model III 


33T = - R p i h q o) - R p i ®J o) 

d Q i 

dT = - Rp i Qi- k P Q i ( p 2 + p 3 + *'* ) “ k P Q i 


dp 


n-1 


n 


+ k 


T 


31T " R P 1 S n-1 + k p ^ p r s n-r ‘ k P p n ^ <0) ' k p p n ^ 


dQ 


n-1 


HF " 

d^ k) 

3 F “ 


£ = k y Q S — k ( R — 1 ) Q p . - k Q m £ o) - k Q M , (o) 
* p Z v n— r p nl pnP pn; 


30=1 

K 

- k p (R - 1) P 1 S Q 0> - k p Sp 1 "' ^ 0)+ k p *1 ^ 


kl ' 

Trrrarr ivi s 


d ^ 


( k ) 


+k P ^ TH&TT 


(l)s(W) 
k 


FF 


- k p <^ 1)P l H Q k> - k p S Q k) " k p fi Q k> ^ 0> + k F £ 

TTT^m ^ k_1) 


dH 


( k ) 


dF— = V^ X) P 


^ k * =( 1 ) V 

i if(k-Tyr *fc + k P f- r 


k» 


1—0 


p £T 0 iTTk^Hi 


nr 


k — 1 


(1) s»(k-l) . 1 „ ^ Jci . 5 (1) 5(k-l) 


+ i k p Z rnrarr *r - y**> 


p i ff s k> - k p fi i k) (^ 0> +^ 0) ) 


Initial conditions? Same as in Table A3.1, 
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CHATTER A-3 

MATHEMATICAL COMPUTATIONS 

In both ARB type and AHA. + ER'B type polycondensations 
we have used the four point Purge Zutta Gill method^ 51 ^ for 
predicting the batch concentrations. Erom these batch concen- 
trations we have calculated the SCSTR outlet concentrations, 
using equations (6) and (8), by Simpson's rule^^, with a 

step size of 0.05 X . For R = 1 .0 when an analytical solution 

(33) 

was possible, a check was made by using laguerre quadrature 
Each of the step by step integration starts at X = 0 and the 
value of each of the integrals (Equations (6) & (8)) is checked 
after every 400 step increments in X. The integration is stopped 
when it is found that even after 400 such steps the maximum change 
in the fastest growing integral is less than 0.59^ of Its present 
value. For ARB type the convergence is decided on the basis the 
integral of the second moment as this is the fastest growing 
integral. For ARA + BR'B type polymerization the convergence 
criterion is based on sum of the integrals for 2nd moments of all 
three species. A check regarding accuracy was made by decreasing 
the step size which did not produce any significant difference. 

After stopping the step by step integration, the actual 

values of the integrals are obtained by extrapolation. This is 

/ 

possible because the rate of growing of the values of each of the 



integrals gradually decreases as ’X’ increases until it becomes 
almost zero. This means that the integral values increase 
asymptotically and the final value is obtainable by extrapolation. 

The integrals for monomer concentrations gives the average 
conversion in the outlet of SCSTR. Similarly average value of the 
first three moment integrals give the number average and weight 
average chain lengths from -which the polydispersity index can 


be calculated. 
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CHA PTES A -4 

RESULTS AND DISCU SSION 

The parameter-' used in the study of ARB type polymeri- 
sation is the ratio of the rate constants , 'R' and in the study 
of ARA + AR'B , are both ’ R' and the initial molar fraction of 

the limiting reactant P, ’ fPo' . The weight distribution of th.- 
chain lengths for ARB polymerisation at -various values of 0 
and the reactivity ratio are shown for models I and II in 
figures 1 and 2. For the first kinetic model, it is observed 
that for R>1, the FWD curves split into two as 'n' decreases: 
one for od'-‘' values of n and the other, having higher weight 
fractions, for even values of n. This is expected because of 
the preferential reaction of monomer with itself thus giving 
a relatively lar.-rr rate of formation of the dimer. The un- 
reacted monomer is very high in this case. In the second 
kinetic scheme, on the other hand, a splitting is obtained 
for R <£ 1 . This can possibly be explained by the consideration 
that the reaction rate constant in dimer-dimer reaction is much 
higher than, say for examnle, between monomer and dimer 

(by which reaction is trimer is formed). These findings are 

(21 22 ) 

same as obtained earlier for batch reactors ’ , and 

HCSTRs^ 27 ’ 28 ) before. 
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The number average chain lengths and polydispersity indices 
for the two models of APB polymerization in SCSTR are plotted in 
fifrs 3 to G as function of 0 . For the sake of comparison, 
corresponding plots for batch and CSTR are included. Also, values 
obtained with fPO = 0.5 for APA + BE'B polymerization in SCSTR 
are included for comparison. The values for the latter ones are 
found to be always less than the corresponding ones in ARB 
polymerization. Usually the curves for SCSTR falls, as expected 
between the curves for batch and HCSTP . 

In figures 7 and 8 are compared the values of polydispersity 

indices for the three different models of ARA +• BR’B polymerization 

when the initial reactant ratio is other than 1:1.’fP0' is the 

original fraction of the limiting reactant. From these plots it is 

seen that (i) the model II is most sensitive to the change in R 

values while model III is the least sensitive; this is shown by 

the fact that the change of is maximum for a change in R values 
to 10.0, 

from R = 0.1. in case of model I. (ii) that in general for such 
a polymerization, the PDI value obtained for R<Cl are much higher 
compared to those obtained for R^>1 irrespective of the model and 
the fPO value. This possibly indicates that according to these 
models, for unequal initial concentrations of monomers, the 
dispersity of the polymer formed will he higher if there is a 
preferential polymer-polymer reaction as compared to the monomer- 
polymer reaction. 
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In figures 9 and 10 cross plots of J* vs monomer conversion 
for ARB type monomers, with constant chain length curves, are 
shown. These curves are of significant interest since they show 
the relationship between the three important parameters, conversion, 
chain length and polydispersity index. Here also it is seen that 
for same conversion f falls as *R' decreases. 

In conclusion it may he observed that though the kinetic 
models that are proposed here are idealized they may he expected 
to simulate the polymerization reaction reasonably well under 
some conditions and for many classes of polymers. The results presen- 
ted here are the outcome of a theoretical investigation and are 
given in the form of dimensionalized plots which can he used for 
study of particular classes of polymerization reactions. At least 
they can he used as a starting step while analyzing for more 
complex aspects of nonequal reactivity, provided more data are 
available. 
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Fig. 7 - Pol ydispcrsity index vs. reactor holding time in 
ARA + BR8 polymerization. 




Fig. % -Polydispersity index vs. reactor holding time in 
ARA ♦ BRB polymerization. 
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CHAPTER B— 1 


Introduction to Nylon-6 

Nylon— 6, which is the generic term for a synthetic linear 
polyamide in which the repeating unit of the chemical structure 
contains six carbon atoms, — ODNH — (CH 2 ) 5 — , is produced by 
the ring— scission polymerization of the cyclic monomer 

caprolactam. The chemistry of Nylon— 6 polymerization has been 
comprehensively reviewed by Reimschueseel * 1 . Polycapro amide 
can be produced from caprolactam by one of three entirely 
different routes viz*, hydrolytic polymerization, catinoic poly- 
merization and anionic polymerization. The cationic and 
anionic cases are only of theoretical interest* but the 

hydrolytic polymerization is by far the most important commercial 
process for Nylon-6 production. The actual plant practice may 
vary depending on local conditions and demands, but the basic 
process remains the same*^ # ^ ^ , Most of the commercial processes 
consist of a feed section where carefully controlled compositions 
are obtained containing molten caprolactam, water, acids or bases 
and other necessary additives. The bulk of the reaction takes 
place in the polymerization reactor, the product from which is 
treated for removal of unreacted monomer and low molecular weight 
cyclic oligomers either by hot water extraction or by vacuum 
evaporation. The monomer recovered is recycled and the polymer 
is sent to the processing section. 

Many variations in the process conditions and the reactor 
designs have been investigated* 1# 3/10 \ Analysis and simulation 
cf different polymerization reactor systems have been carried 
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out by different workers^* and based on these information 

(10 13 14 ) 

some optimization studies have also been made * * • 

Besides the existing processes reported in the 

literature ^ , some novel operations have been patented 

recently ^ 24-29 \ Many process modifications have been suggested 

also to increase the process eff iciency^ 3 '^ \ to improve 

( 1719 ) 

the product properties and for better removal of monomer 

and low cyclics^ 20 " 23 ^ . 

In our present study, all the major established reactions 
involved with the hydrolytic polymerization of caprolactam 
have been taken into consideration . Based on the available 
information on reaction kinetics and polymerization mechanism 
and partly with a semi— empirical approach, equations for three 
sets of molecular weight moments have been developed, correspon- 
ding to the three different species present viz., the bifunctional 
chains, the monofunctional chains and the closed macrorings. 

From these three sets of moments it should be possible to get 
all the necessary information about the monomer conversion, 
molecular weight distribution and their averages for all kinds 
of species at every stage of polymerization. 



CHAPTER B— 2 


Hydrolytic Polymerization of Nylon— 6 

Kinetics and Equilibrium Studies : The course of reaction 

during polymerization of caprolactam can be followed by 

the end group determination as suggested by Willoth (30) * 

Studies have been made on the polymerization kinetics by 

Hermans et.al. (31) ‘,' Kruissink et.al/ 32) ‘, Heikens et.al. <33) , 

Willoth^ ^ 4 ' 33 Dorr and Willoth^ 3& ^, Tai et.al. ^ 37#33 ^ for 

hydrolytic polymerization with or without the presence of 

acids or amines. Studies on polymerization equilibrium 

and effect of process conditions on the MWD of equilibrium 

polymerization have been reported by Fukumoto^ , Majury v , 

40 ( 41 ) 

Tbbolsky and Eisenberg , Reimschuessel , Reimschussel and 
Dege^ 42 ^ and Smith ^ 43 \ 

Mechanism of Polymerization : On the basis of the kinetic and 

equilibrium investigations made so far it has been established 
that the hydrolytic polymerization involves three major 
equilibrium reactions, namely, 

i) ring opening reaction, where the monomer, the caprolactam 
ring, is opened up by water molecules under proper reaction 
conditions, to form amino caproic acid, the shortest 
bifunctional chain, 

ii) poly condensation reaction, where the amino and carboxylic 
acid jid groups in linear bifunctional chains react with 

one another to form larger chains via amide linkages, 
with water as the condensation product; and 



iii) poly addition reaction, where the monomer caprolactam 

adds directly to any of the growing linear chains as in 

44 

a typical chain growth polymerization . Poly addition 
is the primary reaction during the polymerization process 
and is mainly responsible for the conversion of 
caprolactam. But in the later stages of polymerization 
linear poly condensation becomes important since this 
reaction determines the final molecular weight distribution 
in the product, as equilibrium is approached. 

Besides these three, two more reversible reactions 
have been identified. One, the so called termination, is 
the reaction of bifunctional linear chains with monofunctional 
acids or amines, which are most commonly used as chain 
stabilizers^^, to form monofunctional long chains. The other 
reaction, known as cyclization, is a transamidation by which 
the end group of long poly caproamide chains attacks its own 
amide linkages, by what is termed as *back biting' , to form 
shorter linear chains and macro rings of moderate molecular 
weights. This reaction will again be discussed later. 

Cyclic Oligomers ; In the equilibrium polymer of Nylon-6, some 
low molecular weight macrocyclic compounds are always found. 
Though their concentrations are rather low, the removal of these 
compounds from the polymer is essential prior to processing, 
since their presence adversely affects the processing of the 
polymer ard its physical properties. For this reason these 
compounds have been subject to considerable investigations. 



detection 

Their ^ and quantitative estimation has been done by 

( 47 53 'i ( 11 ) 

many workers * . Mochizuki and I to have reported 

cyclic oligomer concentration as a function of reaction time 

during polymerization. Based on the ring- chain equilibrium 

concept of Jaoobson and Stockmayer^ 54 \ Semiyen and Walker^ 55 ^ 

had attempted, using the Flory-Williams rotational isomeric 

state model, to predict theoretically the percentage of low 

cyclic oligcmers in Nylon-6 equilibrium melt, though the 

comparison with experimental data is inconclusive. Andrews 

et al^ 4 ®^ have reported experimentally determined equilibrium 

concentrations of cyclic oligomers up to the hexamer and the 

corresponding equilibrium constants for cyclization reaction. 

They also described the effect of post polymerization treatment 

( 1 ) 

on their concentration. Reimschuessel has attempted to 
curve fit these equilibrium constants as a function of the 
ring lengths and temperature. 

Rate and Equilibrium Constants : Reports on the investigation 

of reaction kinetics of caprolactam polymerization have shown 

that each of the ring opening, poly condensation and polyaddition 

Cl 31 33) 

reactions are catalysed by -CDOH end groups * * , though 

the actual dependence of rates on catalyst concentration is not 
established rationally, in the usual empirical approach the 
forward reaction rate constants of these reactions are thought, 
as a first approximation, to consist of two parts, one of which 
is linearly dependent on -CDOH group concentration, and the other 


independent as follows: 

k± *s k^ + c k| 


(31) 


( 1 ) 



•where 1 i* designates either of the above three reactions* 


*c* is the acid end group concentration; k? and are similar 
to Arrhenious reaction rate constants* each being determined 
by a frequency factor and en activation energy e.g. < 

k^ = A^ exp (- E^ / RT) (2) 

where superscript * j ' is either *o* or ’ c'. The equilibrium 


constants K*s for these three reactions have been determined 


from thermodynamic quantities ' and *As* and many reports 


agree quite well. The relationship is 

As. AH, 

K i = exp t T “ RT ^ 


(3) 


Tb. completely specify the kinetics* one thus needs to know 

A^* E^, A and A.Hj as given in Reimschuessel^ Recently* 

•• 37 55 

however* Tal et al # have determined these values independently, 

using least squares curve-fitting and though their reported 

reported numerical values sometimes differ substantially from 

Reimschuessel 1 s values* they are claimed to produce less 

deviation from their extensive experimental data. 

( 11 ) 

Mochizuki and I to in their model* have tried to 
show that the forward reaction rate constant in the reaction with 
monofunctional acid depends on -CDOH concentration as well as 
the melt viscosity* though their data show that the viscosity 
dependence is quite small. In our case* we will consider that 
the reactivity of an end group in a monofunctional chain* for 
example* is the same as that in a bifunctional chain. 

In the cyclization reaction also* Mochizuki and Ito^ 11 ^ 
have taken an entirely different approach which is not quite 
unquestionable. They have used* for a cyclization reaction* 



the reaction between the -NH 2 and group and -OOOH end group 
of the same chain, with production of one water molecule and 
a macroring of the same length: 

e *3* V==^°n + H 2 0 ' n = 2 ' 3 '*" (4> 

They did not consider the attack on mid chain amide linkages 

by end groups, which is the general cyclization reaction. 

Moreover, their data for the rate constants are not published 

due to proprietory reasons. In the present case we have used 

a different cyclization reaction scheme and as a matter of 

fact, we have not considered the above reaction of Mochizuki and 

I to at all. The reason is if we consider the above reaction 

with the same theoretical approach as in the other reactions, 

it is imperative that we consider the ring opening of 

macrorings also, as it is a reversible reaction. But, except 
(53) 

Tai et.al. who considered ring opening of the dimer, no 

data are available in the literature for such a reaction. 

Hence, we have kept it out of consideration. However, we can 

partially justify it by considering the fact that the 

probability for this reaction would always be the least among all 

other cyclization reactions for the same chain length because it 

involves complete end to end bending of the polymer chain. In 

our model we have been forced to make certain simplifications on 

(54) 

the original model of Jacobson and stockmayer in order to 
make it amenable to mathematical treatment. In a semi-empirical 
approach we have taken the cyclization equilibrium constant to 
vary inversely with the length of the macro ring formed. We 
will consider it in some detail in the next chapter* 



MWD : 


Iheo retical studies on the molecular weight distribution 

( 57 ) 

in polycaprolactam has been relatively fewer. Willoth and 
( SR ) 

Hermans et.al. have tried to match the MWD in poly caprolactam 
to a Flory-Schulz distribution with relatively little success, 
this is more so in the later stages of polymerization. Tirrel 
et.al. simulated the polymerization of caprolactam, in 
different PFR— CSTR combinations, under some simplifying assum- 
ptions (and with some errors) and reported the theoretically 
predicted MWD averages. Tai et.al. recently tried to 
predict MWD averages theoretically with an assumption of 
Schulz— Zimm distribution, and found polydispersity indices 
from GPC about 9% higher than computed valves. 

In our work we will be making use of an approach 
suggested by Krishnan^ 1 ^ and Gupta and Gandhi which is based 

on the method of Min^ 0 ^. Krishnan has simulated some 
CSTR— PFR combinations for Nylon-6 polymerization and has reported 
the MWD averages, while Gupta and Gandhi have suggested a general 
method for getting the actual MWD as well. The mathematical metho 
used by them is computationally efficient, but in none of the 
above works the cyclization reaction has been included, which we 
will be doing here. 

In this context, references to two other recent works 
may be made, though they are not particularly related to Nylon-6. 
In the first, Howell et, al. ' attempts to get a semi empirical 
distribution function for linear and circular Gaussian chains? 
and in the second, Siootte^ uses a Gram-Charlier series 
expansion of molecular weight distribution. This latter approach 



is in many ways similar to that used by us, but comparatively 
more difficult computationally and there is no evidence that 
it is guaranteed to work for Nylon— 6 system* 



CHAPTER 3-3 


Reaction Simulation 


Reaction Scheme : The present reaction scheme which includes 

all five reactions identified so far in Nylon-6 polymerization 
was first proposed by Gupta et.al,^^. The reaction scheme is 
presented in Table B-l. A.t any time in the reaction mixture, 
there would be three kinds of species; e.g . bifunctional chains, 
denoted by n=l, 2,. . . ; monofunctional chains, denoted by 
V n = 1 » 2 '— . and rings, denoted by C n , n = 1,2,.**. 

Initial feed is composed of * C,' (i.e. the monomer); water, 
designated by *W', and 'A^*, the stabilizer (usually monofunctional 
acid, e*g. RCDOH, where * R* is a nonreactive group). The rate 
constants k^, k 2 , k^, k£, k£ and k^ are based on the concen- 
tration of the functional groups like — OOOH, -NH^ and — CDNH- 
etc. and these are the ones reported in the literature^ 1 \ 

In writing rate equations for molecular species, we have to 
account forall the reactions occuring. 

i) for ring opening reaction, rate equations are directly 

obtained for molecular species, e.g, c^* 

... constant 

ii) for condensation reaction, in the forward re ac tion/wou Id 

be 2k 2 , as there are two possibilities, the -GDOH of S n 

may react with -NH 2 of s m and vice versai, unless n and m 

are equal, when it would be In the reverse 

reaction, the rate constant will usually be 2k^, because 

a water molecule can attack a linear chain S„ at two 

n-Hn 



places to give exactly similar products; the rate 
constant will however be if (n-t-m) is even and the 
central anide linkage is attacked* 

iii) for reaction with monofunctional acids, the forward 

reaction rate constant is always k 2 > as only one 

-NH 0 on may react with one -000H of h , The reverse 
z n m 

reaction rate constant will always be k^, because the 

reaction between W and A gives different products 

n+m 

when water attacks the i~th and the (n+m-i)th amide lin- 
kages, 

iv) in the general poly addition step, since the -NHg end 
groups of S n can react with any of the ‘m 1 amide 
linkages in a cyclic m-mer, the rate constant is mk^. 

In the cyclization step, the -NHj ® nc3 group of s n+m 
can attack any of the (n+m— 1) -OONH- groups on itself 
to give several products; the rate constant for this 
reaction depending on the size of the macro-cyclic 
formed (C^) and denoted by k^ (i). 


The cyclization reaction rate constants require 

special consideration here. First of an let us establish 

the relation between this cyclization reaction and that 

( 54 ) 

originally dealt with by Jacobson and Stockmayer , i.e. 


C +- S 
n m 


( 1 ) 


n+m — 

It should be noted that these two reactions are of opposite 
nature; i,e, the way the reactions have been defined, the 
forward reaction of the first is equivalent to the reverse 
reaction in the second, the difference being that the f 
reaction considers equilibrium between all 




A Mil? 



a time, while the latter concentrates on a single chain in 


equilibrium with a ring of a particular size. It has been 
(54) 

shown that the equilibrium constant for reaction (1)', 

K^, defined as 

C nV S r- < 2 > 


K n= 


m 


m+n 


can be expressed as 


K n = V*" 


(3) 


where *x‘ is the fraction of reacted end groups in the chain 
fraction. Now if we consider our present cyclization equilibrium 
it is easy to establish that 


co 
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n+m 
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k 3 (n) 
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c n E 
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(4) 


Defining an equilibrium constant K 3 (n) for the present case 


as K 3 (n) = k 3 /k^(n), we get 


n 


K 3 (n) — n c 


(1 - 


E S m 

m=l 
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“"n 31 ” 


(1 


-ZeSZ 

m=l v _ , 
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) 


m=l 
(5) 


It should be noted that the concentrations and 

•x' are the equilibrium values. Hence x is very nearly unity, 
So from eqn, (5) it may be concluded that at least for small 
values of 'n* , it may be written, wi tiout much error/ as 

1 


K 3 (n) 


n K 


( 6 ) 


n 


Now for n s= 1, is the equilibrium constant for the reaction 
3 n+l* 


+ S n — ^ S_ . This value has been investigated by many 



( 1 ) ' 

workers and agree reasonably well . 

The values of K n for n = 2 to 6 are reported in 
the literature ^* 4 ® \ values for n y 6 are not available because 
it‘s not possible to detect and qu.an ti tatively' measure .the 
concentration of C n in equilibrium polymer beyond n 6.. 

From the data of vs. n, it is shown that fortunately the 
functionality can often been expressed as 

K n = — ^2 / n = 2,3,.... (7) 

.... .. n 

whore B isa constant. As an example *B* is approximately 
0.1692 at 235°c (64) . Thus from (6.) and (7), one gets 

k * B 

k 3 ^ = “^n * n -.2/3/... (8) 

% * 

We have assumed that this functionality between k^(n) and n 
is retained for values of n greater than 6. As will be 
seen later, this functionality is a very convenient form which 
enables us to adopt, the moment approach. For any functionality 
relationship other than that in (8), a moment analysis becomes 
extremely difficult. 

Rate Equations ; The reaction mechanism having been established, 
we may now proceed to write the rate equations for various 

species present in the reaction mass. The rate equations for 

, * ‘ 

a batch reactor or PFR are shown in Table B-2. For a CSTR, 
the derivatives on the LHS of these equations will be replaced 
by the difference in input and output concentrations divided by 
the reactor holding time. Because of the reversibility of 
the reactions, it may be seen that the solution 'for 1 any species. 



s a y s n , will require knowledge of the concentrations of- 
all three species A,c and S, longer than *n*. For this reason, 
the analysis becomes, at least theoretically, the solution of 
an infinite set of equations. Were it possible to solve all 
of these coupled differential equations (or algebraic equations, 
as the case may be), all the necessary information about the 
polymerization could have been obtained. Gupta et.ai.^^ 
have integrated this set of equationSfor batch reactors and 
have reported the results in the form of molecular weight 
distributions and their averages, for both chains and rings, 
as function of reaction time. However, as is expected, this 
approach takes disconcertingly a large computation time, vhtich 
is quite natural because of the necessity of handling of literal 1 
hundreds of coupled differential equations. This has motivated 
search for more computationally efficient mathematical approaches 
even at the expense of some accuracy, if necessaryj The moment 
approach vhich will now be described is undertaken with this 
view in mind. This approach will presumably be more useful, 
in case of CSTR’s, since it will be extrenely difficult, to 
solve so many algebraic equations efficiently. 

Moment equation ; As has just been stated, to make a rigorous 
and accurate analysis of the polymerization reaction, one needs 
to handle, at least theoretically, an infinite number of 
equations. The idea in the moment approach is to convert this 
infinite set of equations to a set of reasonably small number 
of moment equations so that the analysis becomes easier. 



As has been previously described there are three kinds 


of species in the polymerizing mixture; bifunctional chains (S ' 
mono functional chains (A^ ) and rings ( ) ; and each of these 
species will have a weight distribution of its oWri, 

Let us define three sets of moments as follows; 


CO 

- E 

n=l 


n k s n ,1 t = z: n k c n 

n 1 k n 

n=i 


GO 


and 

A 

?> k 

V 


E 


n=l 


n^ A 


n 


(9) 


If we multiply the rate equation for any species, say 
lc* 

s , by n , and sum up the resulting equations for n=l to oo, 
we should get the rate equations for d /\^/dt. Ihe same 
method applies for the A*s and C n ’s, After considerable 
mathematical manipulations one arrives at the set of equ ations 
given in Table B~3, Note that no expression is available 
for d ^ /dt, because that involves summation of a harmonic 

° a*, aft ay 

series. Also note that = 0 ,which is 

expected since ^ ^ « constant = C 10 + a c * 

To integrate the rate equations for the moments it is 
necessary to get explicit rate equations instead of those 
involving series summations as shown in the above equation. 

In order to do it, recurrence formula for these series summation - 
are necessary. Appendix B-l lists some of these series sums. 

Samples of individual moment rate equations are shown in 
Appendix B— 2, Using the series summation formula listed in 
Appendix B-l, any other moment rate equation can be similarly 
obtained. 



It should be noted that the rate equation for the k-th 
moment requires the knowledge of the (k+l)th moment, our 
subsequent approach is to search for a mathematical tool using 
Which one can express the (k+l)th moment in terms of the k-th 
moment, and get the actual molecular weight distribution from 
theee moments. This method will be (described in the next 
chapter. 

Once we are in a position to express the (k+l)th 
moment in terms of the k-th moment, we can integrate the above 
set of differential equations until equilibrium is reached. The 
results thus generated will give all the necessary information; 
and after a sound mathematical method is established, one may 
possibly use this technique for an efficient reactor simulation 
purpose and for subsequent optimization studies. However, 
that is beyond our scope. 



Table B— 1 


Kinetic scheme for nylon 6 polymerisation 


1. Ring opening: 

2 . Polycondensation : 

forward reaction : 

reverse reaction: 


C., + W — =*• S, 

1 -sr 1 


2k-(nj4n) or k_(n=ra) 

s_ + — ^ s + 

n m n+m 


2k) 


S + W 
n+m 


s + s . 

r n+m-r 


r — 1 / 2 , . . . 


Note : The above rate constant will 
be ki if (n+m) is even and 
for r = ^-(n+m) 

mk„ 

3. General polyaddi tion : S n + C m S n+m , n,m « l,2 /# .. 


' Cyclization : 


kl( r ) 

S — ^ S + C , 

n+m n+m-r r 

r = 1,2,... 

Note: The above two reactions are 
also valid between A-species 
and C -species 


4. Reaction with mono- 
functional acid : 


S + A 
n m 


-5* ■ A + W , 
n+m 


Reverse : 


n,m = 1,2,... 

^2 

Vw. + M S r + A n«-r ' 

r = 1,2, . .. 


Contd 



Table B— 1 (continued) 
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Table B-2 


— at ions In isothermal nylon 6 reactor 
^ brackets* ^ Goncentration terms are not enclosed by square 
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dt 
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n-1 cxq 

s n S' MO*) + k 4 <m) S - k S ^ A 

n m=l 3 m=l 3 n+m 2 n m 


dC i 

dt 


dC 

n 

&r 


+ kJ,W 
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\ 
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n = 2, 

3/- 4 


c-C 
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n 

n=i 

+ k' (1 ) 2- S 

3 „ o n 

n=2 

“ k 3 C l *2. 

n=l 


n 
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» 

2 x 



n 
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nC 
* n 

2 s + k* (n) 
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A 


m 
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Table B—2 (continued) 
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Table B~3 


Moment rate equations for nylon 6 polymeria a tlon 
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CHAPTER B-4 


Computational Procedure 

We have seen in the last chapter that the solution for 
the Nth moment requires the knowledge of the (N+l)th moment. 

Our next task is to (i) find the (N+l)th moment from the lower 
moments and (ii) to interpret the N moments in order to get 
the actual molecular weight distribution. We will adopt Min's 
approach^ 0 ^ as has been done before by Krishnan^^ and Gupta 
and Gandhi^ 61 \ 

The starting point in this method is the consideration 
that any continuous normalized distribution G(x ) can be expre- 
ssed by the sum of an infinite number of Laguerre polynomials 

as follows* 

co 

G(x) = ffefj C n L n (x) (1) 

ix 

where (x) is the generalised Laguerre polynomial expressed 

as n . 

£ '*> - H «> 3 Tnr^Wr 4 < 2 > 

j=0 

and A is a constant parameter. 

The reason for using Laguerre polynomials is its more 
general orthogonal properties and the capacity of taking care 
of fractional moments as suggested by Min^ 60 ^. Though theoreti- 
cally an infinite number of polynomials are necessary for ail 
practical purposes* it is possible to express G(x) by a reason- 
ably small number of polynomials. Here we have taken the upper 
limit of summation to be N mgx - 9. The actual choice of N mgx 
depends on the storage capacity of the computer* the accuracy 



desired and the computation speed wanted. Once a finite number 


is used, because of the orthogonal properties, the unknown 

<0( 

coefficients C^, n = 0 to N max , can be computed if the first 


N. 


max inoments hhe normalized distribution G(x) are available 


a 


n 


U&+11 ni / 1 \ J 

TWouTT Z_ ( “ 1} 

j=o 


r (n +K+ 1 ) __ 

r(J 4 *+l)(n-j 7 T 


(3) 


Thus if the first moments of G(x) are known, 

<* 

C^s can be calculated therefrom and the actual 
distribution, G(x), can be expressed with reasonable accuracy. 
Once this is done, the (N max + l)th and highermcments of the 
normalized distribution can be obtained using the general 


expression 


N 


F(pc+ k + 1) 

“rnrm 


max 


n 


+ TTft+rr 21 

n=l 


of 

r-» 

^n 


2 c-i>J 


j=0 


P Cn +c<+ 1 ) P(k+j +^+ 1 ) 

(n - jrm twztt) 


(4) 


for derivation of (3) and (4) see Appendix (B,3) 

The basic approach is thus to convert a given set of 
say *k' original, unnormalized moments to a set of ' k* 
normalized moments. Next by the method described above, the 
higher normalized moments are obtained, from which is 
regenerated the (k+1 )th moment of the original distribution, 

lb convert the original moments to a set of normalized 
moments, the first task is to simulate the actual discrete 
distribution by a continuous distribution. Let us designate such 
a distribution by M(v)dv, which has to be converted to normalized 
distribution G(x) dx. Obviously it follows that 



OD 


M(v) dv = M_i°^G(x) dx 

CO v 


(5) 


since J G(x)dx is unity, 
o 

Min has suggested that there be a linear relation 
between v and x e.g. 

v = a,x. (6) 

the constant 'a' has to be determined from other conditions. 

Since is a free parameter in ©3.(3), it is usually customary 
<X 

to put C x as zero. In such a case it follows from eq.(3) 


that 


G ± = 1 +<X 


Now from eq.(5) we can write 

CD 

M(v) dv = F. 

V 


CD k 

J v M(v) dv = 


CD 

R (o) - / 


or 


Putting k 
a 


a k 3^ 

v v x 

1 in eq. (8)- we get 

r(D 


(7) 


v G(x) dx for any value of 1 k 


( 8 ) 


FT 


M 


( 1 ) 


(!+«) H. 


I oT 


and 


R (k) = f (o) 


V 

r 


V = 




(1+G0 


(l+oO K^ o) 


(k) 

G x 


(9) 


( 10 ) 


Eqn.(10)relates the original moments to the normalized ; 

| 

moments. It should be noted that this eqn. works only for k=2 I 
and above f For ks=l eqn»(10) simply becomes eq. (7). It only 
shows that F00 and G^ are inter-related but neither has a ! 

unique value. To fix a unique value of * o('/ Min has suggested thi 

a .. 

Cj he also put to zero. For such a case, eqn. (3) yields j 


G 2 = (1+00 (2 +ft) = G ± (1 + G x ) 


( 11 ) 



On the other hand putting k = 2 in eq.(10) yields 

fi i 2) 4 0) G 2 

inTT- - ~2" 

M v G 1 

Combining eqs. (11) and (12) one gets 


( 12 ) 


G 1 - 


g(2) ^(o) 

v V 


(13) 


I 


(I) 2 


- 1 


Ihus G ^ and hence get uniquely 1 determined* 

Once the G(x) distribution is obtained, M(v) distribution can 
also be obtained by making use of equations (5) and (9)/ yielding 


M(v) = 


M^ o) (l+C( ) 
“ ^ 


I 


G(X) 


X = 


Hi o) (l+0C) 

— xrr- * v 


H 


(14) 


Thus for integral values of v, M(v) can be calculated; 
and though, strictly speaking, M(v) is a density distribution 

function and not the actual distribution, i,e, actual discrete 

b» b* 

distribution of say S n is related as ^>~ S n = J M(v) dv* 

n-a a 

to calculate individual concentration S n , the evaluation of 
M(v) at that value of v is sufficient. Because Av = 1 in the 
integral, and 4M(v) can be taken to be equal to S n * 
remarks apply to all other distributions* 


Ihe same 



CHAPTER B-5 


Computational Problems and Discussions 

There are two main computational problems present 

in this case. Firstly, it was seen that though the value 

of a suggested by Min gives faster computation, it does 

not necessarily give good convergence. The convergence can 

be best tested in the following way. While using Min's 

procedure we start initially from a set of N max number of 

unnormalized moments from which G-^'s are calculated* Now 

wo may recalculate the normalized moments by putting k = 2,3,... 

...,N in Eq. (4 ) in Chapter B-4, which, when unnormalized, 
max 

can be compared to the original set. The sum of the deviations 
from the original moments constitutes the best convergence 
criterion. The same has been suggested by Sicotte also. 

For a particular set of original moments, this sum of 
deviation is a function of a only, and it was seen that in 
most of the cases, the a suggested by Min gives a poor 
convergence. The value of a that gives the best convergence 
may be obtained by an optimization process, where the sum of 
deviation is to be minimized. This a is subsequently used to 

compute the ( N maX +1)th moment. 

The second problem comes while applying Min's method to 
the ring length distribution. It should be noted that there 
is no way to directly find the zeroth ring moment, since 
it necessitates the summing of a harmonic progression, for 



oo 


which there is no convenient formula. But the value of 
zeroth moment is essential in order to apply Min‘s method. 

Hence one has to think of some way to circumvent these 
problems. 

In the present case we have used <x as a truely free 
parameter for the chain distribution and obtained a by an 
optimization as justified above, one surprising finding is tha' 
integer values of a always give better convergence. This 
practically makes it an integer optimization process. However, 
it is felt that more studies are warranted in this direction. 

In case of rings, a has again been used as afree 
parameter. G^ Is obtained from Eq.(7) and the zeroth moment 
calculated from by Eq. ( 13) , followed by the optimization 
process some as in case of chains. The difference in the 
procedures for chains and rings is that in the former G^ 
is calculated from the zeroth moment and in the latter, it 
is just the opposite way. This is possible presumably 
because there are two inherent degrees of freedom in this 
computational scheme* 

Using the approach outlined above we have tried to 
solve the thirty coupled differential equations which arise 
in the reaction simulation viz. the rate equations for the 
zeroth and first nine moments of both the chains species 
and rate equation for monomer and the rate equations for the 
first to ninth ring moments using a four-point Runge-Kutta-GIll 
method^ 66 ^ with an initial time increment of 0,1 hour. The 



59 


rate and equilibrium constants of Re in schue sse 1 ^ ^ are used 
for an initial water concentration of 0*44 mole/!kg and 
temperature 235°C. 

The molecular weight distribution for the chains can 
bo found by the method outlined in chapter four and they 
match within some errors with the data given by Gupta et al . 

The some method to find the ring length distribution for the 
macrocyclics was not very successful, however, particularly 
at the low values of the polymerisation time. This is 
presumably because the monomer concentration is so much 
more than the total concentration of macrocyclics that the 
distribution is not very well defined (it should be noted that : 
in the calculation of ring moments the monomer concentration 
is also included). This problem may be circumvented by a 
proper scaling procedure in which the monomer concentration 
is deducted from each of the ring moments. This, understandably 
won ' t instill any significant error because the contribution ; 

by the monomer concentration is the same to all of the ring 1 

moments. By this scaling process a reasonable distribution 
may bo obtained for which the concentration for n~L should be i 
ideally zero. In fact, minimization of the concentration of j 
n=l after scaling in itself may be selecting criteria for a , 
but in this case the least squared deviation may be poor. It j 
should be admitted that we are yet to establish a really 
sound selection criterion. 
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Appendix A-1 


Definition 
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At't'BNDIX B-2 


Bifunc tional 
moments 


Monofunction 

Moments 


Sample rate equations for the moments 
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+8 fl^8 + fo^9 ) ’* k 3 ( 1c^ “ B)(8 f7 - 28 p 6 + 56 f 5 “ 70 ^ 4 

+56 (V 28 p2 + 8^ 
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APPENDIX B — 3 


MIN'S METHOD 


G(x) 


a x max 

H“ + i> Jo 


oc a 

L (x) C 

n K n 


where L 


(x) = ^E. (-1)^ P(n+tx + l) 
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a _ x 

x o 
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or rearranging/ 



